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cK"' Abstract 

as 

We developed numerical code for calculation of the extragalactic component 
I of the spectra of leptons, nucleons and ^ -rays resulting from '^top-down" (non- 

Q ' acceleration) models for the case of uniform and isotropic source distribution. 

We explored two different classes of "top-down" scenarios: the wide earlier 
investigated class of X particles coming from collapse and/or annihilation of 
cosmic topological defects (such as cosmic strings, monopoles, etc.) and the 



> 



lO i models of super-heavy long-living X particles with lifetime of the order or much 



o 



greater than the current Universe age. 



0\ . I. INTRODUCTION 



The uhra-high energy cosmic ray (UHECR) events observed above 10^° eV are 
difficult to explain within conventional models involving ffist order Fermi acceleration of 



^ ■ charged particles at astrophysical shocks 0. It is hard to accelerate protons and heavy 

■ nuclei up to such energies even in the most powerful astrophysical objects 0] such as radio 
galaxies and active galactic nuclei. Nucleons of the energies above ~ 4 x lO^^eV are strongly 
decelerated due to photo-pion production on the cosmic microwave background (CMB) — 

^ I the Greisen-Zatsepin-Kuzmin (GZK) effect — on the distances less than ~ 100 Mpc 0. 

■ Ultra high energy electrons lose their energy even faster due to synchrotron radiation in 
galactic as well as extragalactic magnetic field. Heavy nuclei are photodisintegrated in the 
CMB within a few Mpc []^. There are no evident astronomical sources of UHECR within 
~ 100 Mpc of the Earth in the frame of conventional acceleration models. 

To avoid these difficulties one may assume that UHECR are created directly at energies 
comparable to or exceeding the observed ones rather than being accelerated from lower 
energies. In these so-called "top-down" scenarios 7-rays, leptons and nucleons are initially 
produced at ultra-high energies by the decays of supermassive particles generically called X- 
particles. The sources of X-particles could be topological defects (TD) such as cosmic strings 
or magnetic monopoles that could be produced in the early Universe during symmetry- 
breaking phase transitions envisaged in Grand Unified Theories. In the inflationary Universe, 
the relevant TD could be formed in a phase transition at the end of inflation. Another 
possibility is to produce super-heavy long-living particles thermally during reheating epoch 
of the Universe [H] or from vacuum fluctuations during inflation B . 
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There were several papers dedicated to the calculation of the observable spectra of 
UHECR |10|,|T3|,|T^ . The most detailed discussion of propagation of nucleons, photons and 
electrons is given in ref. |T^. It has, however, some drawbacks which we discuss below. We 
have developed an independent numerical code for calculation of the extragalactic compo- 
nent of the spectra of leptons, nucleons and 7-rays resulting from an X-particle decays in the 
top-down models of UHECR, for the case of uniform and isotropic source distribution. Our 
general strategy is similar to that of ref. |Ty], but we use another numerical scheme for the 



solution of transport equations, and take into account process of proton-photon scattering 
Plb — ^ VI disregarded previously. 

In this paper we present the results of numerical calculations of the observable spectra 
in two top-down models. The first one belongs to the class of models examined in |T^ , and 
we use it for testing our code. The second model is based on the decays of superheavy relic 
particles; it was not considered in ref. We shall discuss the differences in the observable 
spectra predicted in different models, and derive constraint on X-particle lifetime in the 
second model. 



II. INJECTION SPECTRA AND EVOLUTION 

The generic top-down mechanism of production of UHECR looks as follows: the X- 
particles with typical GUT masses mx of the order of 10^^ — 10^^ GeV decay into leptons 
and quarks. The strongly interacting quarks fragment into jets of hadrons resulting in 
mesons and baryons. Mesons, in turn, decay into photons, electrons and neutrino. Together 
with baryons these particles form the primary injection spectrum. 

Injection spectrum for a given particle species (a) from the hadronic channel can be 
written as 

(E t) = ^"•'^'^(^) "■g + ^e dNa{x) . 
" ' dt mx dx ' 

where x = E{nq + nf.)/mx, and dNa/dx is the effective fragmentation function describing 
the production of the particles of species (a) by the original quark. The total hadronic frag- 
mentation spectrum dNh/dx in the energy range of our interest is not yet known definitely. 



Some possibilities were described in ref. |T^. In our current calculations we have used the 
following fragmentation function. 



dNhjx) ^ j -xf if xq < X < 1 

dx 1 otherwise ' ^ ' 



where the value xq corresponds to a cut-off energy ~ 1 GeV. Equation (§) is one of the sim- 
plest QCD-motivated forms of hadronic fragmentation spectrum |]T5[. Assuming a nucleon 



content of ~ 3% and the rest equally distributed among the three types of pions, we can 
write the fragmentation spectra as [^0 ,^1 



= (0.03) , (3) 

dx dx 

dN^+ _ dN^- _ dN^o _ /0.97\ dNh{x) 
dx dx dx \ 3 / dx 
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From the pion injection spectra one gets the resulting contribution to the injection spectra 
for 7-rays, electrons and neutrinos. 

Consider now the dependence of the decay rate of X-particles on time. In a wide class 
of models involving topological defects the rate of production (and, therefore, the rate of 
decay) of X-particles can be parameterized as 



dnx_ 4+p 
dt ^ ' 



(4) 



where the index p > depends on the nature of topological defects from which X-particles 
are produced. For example, production of X-particles by a network of ordinary cosmic 
strings in the scaling regime would correspond to p = 1 if one assumes that a constant 
fraction of the total energy of closed loops goes into X-particles [pT] , p!2| . Models based 
on annihilation of magnetic monopoles and antimonopoles |1T5|JT6| predict p = 1 in the 
matter dominated and p = 3/2 in the radiation dominated era whereas simplest models 
involving superconducting cosmic strings lead to p = ||18|. A constant comoving injection 
rate corresponds to p = 2 and p = 5/2 during the matter and radiation dominated era, 
respectively. 

Another model that we will consider involves primordial X-particles with the mass of 
order 10^^ eV and the lifetime of the order of the age of the Universe or greater. These 
particles may be produced thermally during reheating epoch p| or from vacuum fluctuations 
during inflation [^] and may constitute a considerable fraction of Cold Dark Matter (CDM). 
X particle density in this model is 



nxit) 



[l + zf -i^ 



rii 



[I + Zi) 



(5) 



where rij is X-particle density at some arbitrary initial moment ti^Xi z is the redshift corre- 
sponding to time t, Zi^x is the initial redshift, and lO^^years < tx < lO^^years is X-particle 
lifetime. We normalize X-particle concentration by assuming that at z-i ~ 100 — 1000 (matter 
dominated phase) the density was approximately equal to the critical density at that time. 
Then we get 



Hi 



nix 



'l + Zi 



(6) 



where po is today's critical density. Taking into account that ti <^t (the major contribution 
to the observable flux comes from z < 10) we obtain 



nx{t) 

Finally, X particle decay rate reads 



mx 



[l + zfe ^x. 



dn, 



'•^ = (1 + zf 



dt 




nx{t) 
l + z)- 



Po 



nixTx 



zfe . 



(7) 



(8) 



For the full determination of the injection spectra one needs also to specify the X-particle 
decay mode. In what follows we assume that X-particle decays on nq quarks and rig electrons 
and all these initial decay products share energy equally: Eg = Ee = mx/iriq + Ue). 
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One remark concerning the model of primordial X-particles is in order. One has to take 
into account the fact that these particles are concentrated mostly in galaxies, and so the 
electron component of injection spectra is heavily suppressed due to synchrotron radiation in 
the galactic magnetic field. For the same reason there is an additional contribution in photon 
spectrum. In the case of the decay mode X — > Uqq typical electron energy is 10^^ — 10^'' eV 
for TUx = 10^^ eV and Ug = 2, 3. In the galactic magnetic field of strength B ~ 10~^G the 
energy of the synchrotron photons produced by such electrons is about 10^^ — lO^^eV. Photon 
attenuation length at those energies is extremely small due to pair production on microwave 
background, so one should expect that the synchrotron component is also negligible. In the 
present paper we approximate this effect by simply excluding electrons from the injection 
spectra. It is worth noting that the contribution of our Galaxy to the observable flux may 
be non-negligible. In this paper we present calculation of the extragalactic component of the 
observable flux only. The calculation of the Galactic contribution requires separate detailed 
calculation |[T9[] . 



III. INTERACTIONS OF ELECTRONS, PHOTONS AND NUCLEONS WITH 
BACKGROUND RADIATION AND EXTRAGALACTIC MAGNETIC FIELD 



The 7-rays and electrons produced by X particle decay initiate electromagnetic (EM) cas- 
cades on low energy radiation fields such as the cosmic microwave background (CMB). In 
the absence of magnetic field EM cascades are driven mostly by the cycle of pair production 
(PP; e~e^) and inverse Compton scattering (ICS; 67^ — »■ e'7). The energy degrada- 

tion of the "leading" particle in this cycle is slow. Other EM interactions that influence the 
7-ray spectrum in the energy range of interest are triplet pair production (e7{, — ee~e+) 
and double pair production (77^ — > e~e^e~e"'"), as well as synchrotron losses of electrons in 
the large scale extragalactic magnetic field (EGMF). 

The relevant nucleon interactions implemented in the code are pair production by protons 
{pib pe~e^), photoproduction of single or multiple pions (A^7fe — mc, n > 1) and 
neutron decay. Production of secondary 7-rays, electrons, and neutrinos by pion decay is in 
general negligible in the scenarios where injection is dominated by 7-rays and leptons over 
nucleons. Apart from these processes described in detail in ref. [|10|, for protons of energies 
10^^ — 10^^ eV the inverse Compton scattering must be taken into account. Its cross section 
is identical to that for electrons (with the replacement of me by m^). 

There are two major uncertainties in the parameters which influence the particle trans- 
port. The first one concerns the intensity and spectrum of URB for which there exists only 



an estimate above a few MHz frequency p2[. We adopt the so-called "minimal" URB es- 



timate that has a low-frequency cutoff at 2 MHz [^. The other uncertainty is the mean 
value of EGMF and its time evolution. Conventionally, concerning particle propagation, 
EGMF is assumed to be constant and lying in the range lO^^^G < B < IQ-^G 
The other possibility is to assume that the total energy of EGMF is constant; this implies 

5 -(1 + ^)3/2. 

We assume a fiat Universe with no cosmological constant, and a Hubble constant of 
Hq = 0.75 in units of 100 km sec~^Mpc~^ throughout. 
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IV. SOLUTION OF THE TRANSPORT EQUATIONS 



An example of the transport equation for electrons which includes for simplicity only 
pair production (PP) and inverse Compton scattering (ICS) looks as follows: 

j^N,{E,,t) = -N,{E,,t) J den{e) J dfi ^ ajcsiE^e, fi) + (9) 

/ dE'ME',,t) I denie) J df^^-^^i^^'^ e, f^) + 
J dE^N^{E^,t) J den{e) J dfi ^ ~^ ^^^^ {E^; E^, e, /i) + Q{Ee,t), 

where Ne{Ef,,t) is the (differential) number density of electrons with energy E^ at time t, 
n(e) is the number density of background photons of energy e, Q{Ee, t) is an external source 
term for electrons of energy E^ at time t, fi is the cosine of the interaction angle between 
the CR electron and the background photon {fi = —1 for a head-on collision), and jSe is the 
velocity of the CR electron. The first term in (^ describes the loss of electrons due to ICS, 
the second one accounts for influx of electrons scattered into corresponding energy range due 
to ICS, and the third term describes the influx of electrons produced due to PP by photons. 
The processes with continuous energy loss dE/dt = —f{E, t), such as synchrotron radiation, 
may be accounted for by including into the transport equations the additional term 

dN ,dN , , 

In order to solve the transport equations numerically, we divide each decade of energy 
into n ~ 20 equidistant logarithmic bins and define the central value of the i-th bin as E^. 
Then we replace the continuous integrals by finite sums, and integrate Eq. (|]) over one CR 
energy bin. The resulting set of differential equations has the form 

= -K< + E E PiiNi + Ql- (11) 

where subscripts a and c refer to particle species. The coefficients a^, and depend 
on time and may be easily obtained from the original transport equation. For numerical 
solution of (|TT]) we used the implicit scheme with Richardson extrapolation [ P^ . 

To account for redshifting S.Lee in ref. [|l^] performed the operation Na{E,z) — >■ [1 + 



Az/{1 + z)]~^Na {E[l + Az/{1 + z)], z) for each particle species a after a step Az in redshift. 
Here Az was matched conveniently with the logarithmic energy bin size, log^o[l + ^-^/(l + 
z)] = logiQ^Ei/ Ei^i) = 1/n, which corresponds to the transformation N^{z) — * [1 + A2;/(l + 

zr'K-^'iz). 

It is important to note that the time interval At corresponding to Az may be rather 
long. For example at 2 ~ and n = 20 one has Az = 0.12, which corresponds to distances 
Al ~ 360//ioMpc ~ 500 Mpc and is much larger than the synchrotron loss length of electrons, 
D ~ IMpc, at ~ lO^^eV in the magnetic field B = lO^^^G. So, we have to divide time 
interval into smaller steps to account for all the fast processes such as synchrotron losses of 
electrons, GZK-cutoff for nucleons, etc. This was taken into account in the code. 
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V. RESULTS 



In figs. |T]-^ we present the results of our flux calculations. Since we have not taken 
into account the propagation through our Galaxy, the spectra presented refer to the fluxes 
of particles near the boundary of the Galaxy. Even in the absence of galactic sources, 
the observable flux of electrons will be smaller by several orders of magnitude at energies 
E > lO^^eV, while the photon flux will contain additional component due to the synchrotron 
radiation in the galactic magnetic field ||24|| . 

We simulated the situation of fig. 14(a) in ref. |^ (see fig. §). As a result we got the 
spectrum whose shape is mostly similar to fig. 14(a) in |^5||. But there are some discrepancies 
in details which concerns especially the proton component and may in part be attributed 
to including elastic proton-photon scattering and also slightly different accounting for pair 
production by protons. 

The fluxes of UHECRs in the model of long-living X-particles are proportional to r~^, 
which allows to impose constraints on X particle lifetime. For example, in the case of decay 
mode X ^ qq and mx = lO^^eV, comparison of the levels of predicted spectra with the 
observed one leads to the constraint 

7^ > 10^%, (12) 

where to is the today's age of the Universe and fix is the relative density of X-particles at 
z = 100 - 1000. 

The main difference in spectra resulting from TD and primordial X-particle models 
concerns the ratio of nucleon and electron-photon components at energies above lO^^eV. 
The reason is the absence of electrons in the effective injection spectra in the second model. 
We expect that this difference in spectra will persist even after including into consideration 
the galactic component of the UHECR. 
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FIG. 1. Extra-galactic component of the spectra in the model of long- living X-particles with 
r = lO^to (^0 is the age of Universe) and the decay mode X — > qq. The strength of EGMF is 
B = 10^^^ G. Synchrotron losses in the Galaxy are not taken into account (see comments in Sect. 
0). 
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FIG. 2. Same as fig. || but for B = lO-^^G. 
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FIG. 3. UHECR spectra near the boundary of our Galaxy in the top-down model with p = 1, 
Mx = lO^^eV and constant EGMF B = IQ-^^G, assuming decay mode X — > qq. The spectra are 
normalized to the observed flux at energy 
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FIG. 4. Same as fig. |3| but for the decay mode X qe. 
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